o 
o 

< 

(N 



CD 



c^ 



O 
> 

in 

m 

o 
\o 
o 

-(— > 



I 

o 
o 



X 



Measuring fractional charge and statistics in fractional quantum Hall fluids through 

noise experiments 

Eim-Ah Kim,^'^ Michael J. Lawler,-'^ Smitha Vishveshwara,^ and Eduardo Fradkin^ 

'Department of Physics, University of Illinois at Urhana- Champaign, 

1110 W. Green St. , Urbana, IL 61801-3080, USA 

^Stanford Institute for Theoretical Physics and Department of Physics, Stanford University, Stanford, CA 94305, USA 

(Dated: February 6, 2008) 

A central long standing prediction of the theory of fractional quantum Hall (FQH) states that it 
is a topological fluid whose elementary excitations are vortices with fractional charge and fractional 
statistics. Yet, the unambiguous experimental detection of this fundamental property, that the 
vortices have fractional statistics, has remained an open challenge. Here we propose a three-terminal 
"T-junction" as an experimental setup for the direct and independent measurement of the fractional 
charge and statistics of fractional quantum Hall quasiparticles via cross current noise measurements. 
We present a non- equilibrium calculation of the quantum noise in the T-junction setup for FQH Jain 
states. We show that the cross current correlation (noise) can be written in a simple form, a sum 
of two terms, which reflects the braiding properties of the quasiparticles: the statistics dependence 
captured in a factor of cos in one of two contributions. Through analyzing these two contributions 
for different parameter ranges that are experimentally relevant, we demonstrate that the noise at 
finite temperature reveals signatures of generalized exclusion principles, fractional exchange statistics 
and fractional charge. We also predict that the vortices of Laughlin states exhibit a "bunching" 
effect, while higher states in the Jain sequences exhibit an "anti-bunching" effect. 



I. INTRODUCTION 

Bosc-Einstcin statistics of photons and Fermi-Dirac 
statistics of electrons hold keys to two major triumphs 
of quantum mechanics: the explanation of the blackbody 
radiation (that started quantum mechanics) and the peri- 
odic table. The spin-statistics theorem states that parti- 
cles with integer(half-integer) spins are bosons (fermions) 
and that the corresponding second-quantized fields obey 
canonical equal time commutation (anticommutation) re- 
lations. In three spatial dimensions, the spin can only be 
integer or half-integer since the fields should transform 
like an irreducible representation of the Lorentz group 
50(3,1) (rclativistic) or 5*0(3) (non-relativistic). Conse- 
quently particles in three dimensional space(3D) have ei- 
ther bosonic or fermionic statstics. In one spatial dimen- 
sion (ID) on the other hand, neither fermions nor hard 
core bosons can experience their statistics since they can- 
not go past each other. As a result, statistics is essentially 
arbitrary in one spatial dimension (which in some sense 
can even be regarded as a matter of definition.) In par- 
ticular, the excitations of (integrable) one-dimensional 
systems are topological solitons which have a two-body S- 
matrix which acquires a phase factor upon the exchange 
of the positions of the solitons. In this sense, one can as- 
sign an intermediate statistics to the solitons. However in 
two spatial dimensions the situation is quite different. It 
has long been knownii^ that in two dimensions an inter- 
mediate form of statistics, /raciiona/ statistics is possible. 
A specific quantum mechanical construction of a particle 
with fractional statistics, proposed and dubbed an any on 
by Wilczek— , consists of a particle of charge q bound to a 
solenoid with flux 0, where qcj) bears a non-integer ratio 
to the fundamental flux quantum (J)q. In this paper, we 
study the 2D setting of the quantum Hall system as an 



arena for displaying such fractional statistics, and pro- 
pose a concrete experiment for measuring its effects. 

Not long after the first observation of the Fractional 
Quantum Hall (FQH) effect^, Laughlin proposed that 
quasiparticle(qp)s/ quasihole(qh)s of these incompress- 
ible fluids carry a fractional charge e* — ±i>e deter- 
mined by the precisely quantized Hall conductance, i^ = 
l/(2n + 1) for an integer v^ (for Laughlin states). Soon 
after, it was shown that these qp's, would have fractional 
braiding statistics^, i.e. the two qp joint wave function 
would gain a complex valued phase factor interpolating 
between 1 (boson) and —1 (fermion) upon exchange: 



*(ri,r2) = e^«*(r2,ri). 



(1.1) 



with < 6 < TT. In the anyonic picture, the statistical 
angle corresponds to 6 = ■nq(j)/(j)Q. Laughlin qp's are a 
specific example with q = v. In order to compute the 
statistical angle 9, Arovas et al.— first considered the the 
adiabatic process of one qp encircling another. By calcu- 
lating the Berry phase associated with this process they 
showed that the two qp wave function \I'(ri,r2) gains 
an extra "statistical phase" of e^'^'^'^ in addition to the 
magnetic fiux induced Aharonov-Bohm phase upon encir- 
cling. An adiabatic exchange of two qp's can be achieved 
by moving one qp only half its way around the other and 
shifting both of them rigidly to end up in interchanged 
initial positions. They thus argued that the phase factor 
picked up by a two qp joint wave function upon this ex- 
change process should be precisely half of what it is for 
a round trip and hence 6 = vir. 

From a more general standpoint, anyons are excita- 
tions which carry a representation of the braid group (see 
below). This is possible both for non-relativistic par- 
ticles in high magnetic fields (as we will be interested 
in here) as well as some field-theoretic relativistic mod- 




FIG. 1: An adiabatic exchange between two qp's braids 
worlds lines representing histories of each qp. 



els. As in the FQH example, the quantum mechanical 
amplitudes for processes involving two or more of these 
particles (regarded as low-energy excitations from some 
more complicated system) are represented by world lines 
which never cross (as if they had a hard-core repulsion) . 
In 2 + 1 dimensions different histories of these particles 
can be classified according to the topological invariants of 
the knots formed by their world lines (see Fig.^J. These 
topological invariants are representations of a group, the 
braid group. The statistical angle is one such label and 
it corresponds to a an Abelian (one-dimensional) repre- 
sentation of the braid group. Most quantum Hall states 
are Abelian (in this sense). However a number of non- 
Abelian states have been constructed {e.g. v = 5/2, 12/5 
and a few others). 

Fractional (or braid) statistics is a generalization of the 
concept that fcrmions have wave functions which are odd 
upon exchange (i.e. thus they obey the Pauli Exclusion 
Principle) while bosons have wave functions which are 
odd upon exchange. The notion of statistics is (as its 
name indicates) also related to the problem of count- 
ing states. Some years ago, Haldanei introduced the 
concept of exclusion statistics^ which is a generalization 
of the Pauli exclusion principle for fcrmions. Exclusion 
statistics in finite systems is determined by how much an 
addition of a particle diminishes the number of available 
states for yet another addition. For fermions, the number 
of states will diminish by one if a particle is added to the 
system while for bosons, the number of states would stay 
the same; Haldane's idea was to consider more general 
possibilities interpolating between fermions and bosons. 
Although the Hilbert space counting definition of statis- 
tics and the braiding definition of statistics coincide with 
one another for Laughlin states, the two definitions are 
not equivalent in general. It was argued by Haldane^, 
and shown explicitly by Van Elburg and Schoutena^, that 
the qp's of the FQH states also obey a form of exclusion 
statistics. 

Fractional charge and statistics are fundamental prop- 
erties of qp's emerging from the strongly interacting FQH 
liquid. Each of the filling factors ly displaying precise 
quantization of fractional Hall conductance a^y = ve^ jh 
represent a distinct phase characterized by nontrivial in- 
ternal order- topological order,- which is robust against 



arbitrary perturbations2ii£. The measurement of frac- 
tional statistics will prove the existence of such topo- 
logical order. While there is substantial experimen- 
tal evidence for the fractional charge of these qp's and 
qjjigiij2ji3ji4^ in particular from two-terminal noise ex- 
periments, similar evidence for statistical properties or 
the connection between two distinct ways of determining 
statistics is still lacking. The main challenge in detec- 
tion of anyons is to manipulate these "emergent par- 
ticles" which cannot be taken outside the 2D system 
and to measure the effect of statistical phase. In this 
paper, we show that the noise in current fluctuations 
in a three-terminal geometry (a "T-junction") in FQH 
Jain statesi^ behaves as in a Hanbury-Brown and Twiss 
(HBT) interferometer— with clear and independent sig- 
natures of fractional charge, fractional statistics and ex- 
clusion statistics. 

Our interest in Jain states is two fold. First, the most 
prominent FQH effects which are observed lie in the Jain 
sequence with v = „ ^ . where n and p are integers. 
Further, as shown in the table ^ Jain states include 
"non-Laughlin states" with fractional charge and statis- 
tical angle that are different from the filling factor, i.e. 
e*/e 7^ O/i: ^ v. Therefore one can pursue signatures 
of these fractional quantum numbers independently and 
predict how each fractional quantum numbers will man- 
ifest itself in different aspects of experimental data. Par- 
ticularly, by comparing two FQH states whose qps have 
same fractional charge but different statistics one can fo- 
cus on the effects of fractional statistics. 
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TABLE I: Fractional charge e* and statistical angle 9 for 
qp's of FQH states at filling factor v = p/(2np+l) (n, p 
are inte gers ) in comparison with bosons and fermions follow- 
ing Ref. llTll . Alternative descriptions predict somewhat dif- 
ferent spectrum of qp's for non-Laughlin states at p 7^ 1 (see 
Ref.[T3) and implications of such differences will be discussed 
in the later part of this paper. 

In this paper we present a theory of noise cross cor- 
relations of the currents in a three-terminal junction (or 
T-junction) which we introduced in Ref.|l8l|. where we 
presented a summary of results of the theory. Here, we 
describe the theory and its conceptual and technical un- 
derpinnings in a more detailed and self-contained fashion. 
The bulk of the paper consists of the description of the 
cross current correlation between tunneling currents par- 
titioned from one edge of the proposed T-junction setup 
into two others. The lowest order non- vanishing con- 
tribution to the correlation is fourth order in tunneling 
and it is the quantity which contains signatures of both 
fractional charge and fractional statistics. In particular, 



by analyzing the components of the correlation in depth, 
we pinpoint the role of statistics in various tunneling pro- 
cesses and discuss implications of the connection between 
braiding statistics and exclusion statstics. 

The paper is organized as follows. In Section ^ we 
briefly review the effective theory for edge states devel- 
oped in Ref. [i3| which we will use throughout. In Sec- 
tion lllll we introduce the T-junction setup we arc propos- 
ing and we define the normalized cross current noise to 
be calculated purturbatively in Section llVl In Sect ions Ivl 
and I VII we present the time dependent noise and its fre- 
quency spectrum respectively. We summarized our re- 
sults and discuss its implications for experiments in Sec- 
tion IVIII In four appendices we present a summary of 
the theory of the edge states that we use here (Appendix 
IXj) . of the Schwinger-Keldysh technique and conventions 
(Appendix 0). some useful identities for vertex operators 
(Appendix Q . and the unitary Klein factors (Appendix 
p. 



II. AN EFFECTIVE THEORY FOR EDGE 
STATES 

In this section, we summarize the properties of edge 
states for the Jain sequence and their associated quasi- 
particles. While the creation of qp excitations in the 2D 
bulk of FQH systems has an associated energy gap, {i.e. 
FQH liquid is incompressible), the one dimensional (ID) 
boundary defined by a confining potential can support 
gapless excitationsifliiS. Further there is a striking one- 
to-one correspondence between qp states in the bulk and 
at the edge making the FQH liquid holographioiS. Given 
that edge states comprise an effective ID system prop- 
agating only along the direction dictated by the mag- 
netic field (a chiral Luttinger liquidiStiS) , they can be 
described in terms of chiral bosons using the standard 
bosonization approach. The edge effective Lagrangian 
density for Jain states, derived from the boundary term 
of the fermionic Chern-Simons thcorySS by Lopez and 
Fradkin-7 is given by 
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(2.1) 



where (j)c is a right moving charge mode whose speed is 
set to be Wc = 1 (</)- with g = 1/v and v = \ case of 
Appendix ^1 and ^^v is the non-propagating charge neu- 
tral topological mode obtained as a ujv — * 0+ limit of a 
rightmover with g^q = —1. Note that for the purpose of 
rcgularization and to carefully keep track of short time 
behavior which is crucial for ensuring the correct statis- 
tics, we shall keep a neutral mode speed vm and take the 
limit Vm — ^ only at the very end. 

The normal ordered vertex operator creating a quasi- 
particle excitation on this edge, with fractional quantum 
numbers tabulated in the Table Q] is given bjiii 



V' 



i(T;'t>c+^J-i-+^<pN) . ^ 



(2.2) 



where we define a short hand notation if to represent the 
appropriate linear combination of the charge mode and 
the topological mode 



ip = {-4>c + \ I + -4>n) 



(2.3) 



By noting the equal time commutator for the new field 
if being 



\(p(x, <), fix' , t)] = ITT { -^ 1 I sgn(a: — x') 

\P P J 

= —i9sgn(x — x') 



(2.4) 



the quantum numbers of -0^ can be verified as follows. 
The charge density operator jo = ■^dx4>c measures the 
charge of the qp operator since \jo{x), ■(/'^(x')] = —5{x — 
x')il:i'^{x'). By expanding the vertex operator Eq. H2.2|l . 
this commutator can be calculated order by order to give 



[jo(x),^t(3;')]^^[a^^^(^)^e-f0c(.')] 



1 i 
2t: p 
1 vn 
27r p 



da:[(l)c{x),(l)c{x')] 

dxsgn{x — x') + ■ 



(2.5) 



^ -6{x - x')i}\x'), 
P 

hence e*/e = -. As for statistics, using the Campbcll- 
Bakcr-Hausdorff (CBH) formula for exponential of opera- 
tors A, e^e^ = e^e^e^^'^\ the statistical angle of the qp 
operator can be calculated via exchange of qp positions 
at equal time as the following 



'il)\x)i}\x') = V'^(x')V'^(a;)e-['^(^)''^(^')l 



(2.6) 



which confirms 6* = (1 — —v) modulo 27r. 

Combining charge mode and neutral mode propagator, 
one can obtain the propagator for the field ip: 



{p{x,tMQ,0)) 



^(0c(a;,t)0c(O,O))+ 1 + - (0jv(x,t)(^jv(O,O)) 
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sgn(u7vt - x). 



(2.7) 



From Eq. (|2.7|l . one arrives at the following form for the 
qp propagator at zero temperature (see Appendix [nj: 



(^(x,t)Vt(0,0)) 



= lim 



To 



To + i{t- x) 



(2., 



g-if (l+i) sgn(Djvt-a;) 



r 



Eq. H2.8|l shows that the qp operators have the scahng 
dimension -^ = t^ = t^ztttzt-ttt • Further, as r ^ 0+, 



2 2j?^ 2p(2np+l, 

the hmit a; — > 0^ of Eq. H2.8|l becomes 



(^(0-,i)V'^(0,0)) = 



-K 



^ie sgn(t) 



(2.9) 



with the exphcit dependence on the statistical angle 
emerging in the long wavelength limit. This equal po- 
sition propagator, which explicitly encodes information 
on statistics, plays a prominent role in subsequent point- 
contact tunneling calculations. Most importantly, the 
amplitude of the propagator Eq. (|2.9|) is governed by the 
scaling dimension but the phase is solely determined by 
the statistical angle. 

Since edge states are gapless excitations, they act as 
windows for experimental probes to observe features of 
the FQH droplet. In particular, tunneling between edges 
allows a viable access to single qp properties, provided 
the tunneling path lie inside the FQH liquid (qp's only 
exist within the liquid)^.In fact, such tunneling has been 
successfully used in a geometry hosting a single point con- 
tact and two edges for the detection of fractional charge 
through shot noise measurementsiSiiiii. In principle, 
two edges would suffice for anyonic exchange as well with 
qp's from each edge exchanging positions via tunneling 
through the 2D FQH bulk. However, in order to detect 
the phase information, the path of exchange needs to be 
intercepted, for instance, by inserting a third edge and al- 
lowing qp exchange between first two edges to go through 
the third edge. This third edge thus acts at once as a 
stage and an "observation deck" . Further, if the first two 
edges are each adding/removing qp's to the third edge, 
one should be able to see the generalized exclusion prin- 
ciple in action by probing the third edge. In the next 
section, we propose a "T-junction" interferometer as a 
minimal realization of such a situation. 



III. THE T-JUNCTION SETUP 

In our proposed T-junction interferometer of the type 
shown in Fig|21 top gates can be used to define and bring 
together three edge states I = 0, 1, 2 which arc separated 
from each other by ohmic contacts. Two edges I = 1 and 
I ~ 2 each separately form tunnel junctions with edge 
and upon setting the edge at relative voltage V to the 
two others, qp's are driven to tunnel between edge and 
edges 1,2. 



Denoting charge and neutral modes for edge / = 0, 1, 2 
by cj)''^{xi,t) and 4'^^{xi,t) respectively, each edge state 
I in the absence of tunneling (hence the superscript 0) 
can be described by Lagrangian densities of the form in 
Eq. ITTJ: 



,(0) 



^TTV 



d,,cbi{^dtcf>[ - d^,cbi) + — (a,>^5,0'^.) 



An 



(3.1) 

where xi now represents the curvilinear abscissa along 
each edge, and the total Lagrangian density becomes 



£(0)= J2 A 



(0 



(3.2) 



/=0,1.2 



Boson fields 0[. and (p'-j^ each have the same properties 
(such as commutation relations and propegators) as (f)c 
and (j)N of a single edge described in the Sec. |n] Fields 
of different edges commute with each other, i.e. 



^c/N ' ^c/N 



] = 0, for / 7^ m 



(3.3) 



for l,m = 0, 1,2. Within each edge I, a vertex operator 
e"^' , where (pi again is a short hand notation (p defined 
in Eq. 12.3|l for each edge I, describes a qp excitation 
^j in that edge. However, since tpi's are constructed to 
commute between different edges, we need to introduce 



V2 




FIG. 2: The proposed T-junction set up. Solid (blue) lines 
represent edge states where the edge state is held at poten- 
tial V relative to edges 1 and 2, i.e. Vb — V\ = Vb — V2 = V^. 
Dashed (red) lines show paths of two quasiparticle tunnel- 
ing through a FQH liquid from edge state to edge 1 at 
xo = —a/2 = Xi , xi = 0; from edge to edge 2 at 
xo = a/2 = X2, X2 = 0. The direct tunneling between edge 1 
and 2 is turned off by setting the two edges at equal potential. 



unitary Klein factors Fi with the appropriate algebra, as 
shown in the Appendix^ to give the correct qp statistics 
upon exchange between different edges. Hence, we have 



V' 



\/2naQ 



-.Fi 



^e^^W^+p"^") 



-\/27rao 



Fie'^\ (3.4) 



where ag is the short distance cutoff and the Klein factors 
obey the statistical rules 



FiF„ 



' FrnFl 



(3.5) 



where aim = ~«m; and ao2 = 021 ~ ^c^w = 0- No- 
tice the Klein factors do not affect the statistics between 
operators on the same edge since au = from the re- 
quirement aim = -a„ii- 

Now, if we set the the distance between two tunneling 
points on edge to be a small positive real number a 
(which we later send to zero) and the tunneling point on 
edge 1,2 to be the origin for the abscissa coordinate of 
edge 1, 2, then qp's tunnel at two locations , a:o = — a/2 = 
Xi and x'o = a/2 = X2 on edge and at xi —Q, X2 = 
on edges 1, 2. The operator Vj{t) which tunnels one qp 
from the edge to edges j '=: 1, 2 at time t can then be 
written in the following form: 



V]{t) = MXj.t)i;]{Q,t) 



(3.6) 



The additional term in the Lagrangian due to the pres- 
ence of the tunneling can be written as Lint = imt,i + 
Lj„( jS^iSi with tunneling term involving edges j = 1,2 
given by 



L 



int , J 



{t) = -T,VJ{t) + h., 



(3.7) 



where Tj is the tunneling amplitude between edges and 
j = 1,2. The non-equilibrium situation of setting a con- 
stant DC voltage bias V between edges and j corre- 
sponds to substituting Tj by Fje"*""* in Eq. (|3.7|l . where 
the Josephson frequency is explicitly determined by the 
fractional charge to be ujq = e*V/H (the Peierls substitu- 
^jpjj22j24j25-j Using the Heisenberg equations of motion 
for the total charge Qi = J dxijoi [xi ) of edge I where 
joi{xi) = ^dxcj^i is the charge density operator for edge 
I, as was done by Chamon et al.^'— for Laughlin states, 
one can show that the tunneling current operator, Ij (t) , 
from edge to edge j is given by 



^.W-*^r,(e-"*v;^ 



-iujQt 



v; 



E^r, 
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y. 



(e) 



(3.8) 



e=± 



where we assumed V j to be real and introduced the no- 
tation V'l = !>/+' and V,- = v\~\ i.e., 

V\^\t) = [FoF-^y f.^^Vo{X,,t) ^-^evAO■t)^ (39) 



to represent the summation over hcrmitian conjugates in 
a compact form. Notice that e = + corresponds to a qp 
tunneling from the edge j to and e = — corresponds to 
a reverse direction tunneling. 

We treat the non-equilibrium situation within the 
Schwinger-Keldysh formalis m^^i^?i^^i^^'^i^^i^Ti^? for com- 
puting expectation values of operators. As a simple ex- 
ample, the expectation value of the current on edge j 
takes the following form up to a normalization factor 



{!,) = i ^(TK/,(i'')e'/-^--(*^)'^*Oo 



(3.10) 



n=± 



where 77 labels the branch of the Keldysh contour: ry = + 
for the forward branch (t = —00 -^ 00) and rj = — for 
the backward branch (t = oo ^ —00). Tk indicates that 
the operators inside brackets should be contour ordered 
before taking the expectation value, and ()o indicates the 
expectation value with respect to the free action 5*0 = 
jj^dt J d^xCo where the free Lagrangian density £0 was 
given in Eq. (|3.2|l . 

Our goal is to identify a measurable quantity that dis- 
tinctly captures the role of the statistical angle 9 which 
features in the single qp operator shown in Eq. H2.6|l . The 
normalized cross correlation (noise) S(t) between tunnel- 
ing current fluctuations (A/; ^ h — {h))'- 



{Ah{t)Ah{t')) 



(3.11) 



turns out to be the simplest quantity that can exhibit the 
subtle signatures of statistics and distinguish them from 
those of fractional charge. In fact, in contrast to shot 
noise, this cross current correlation S{t — t') carries sta- 
tistical information even at the lowest non-trivial order, 
which we calculate perturbativcly in the next section. 



IV. PERTURBATIVE CALCULATION 

In this section, we perform a detailed analysis of 
the normalized cross correlation at fourth (lowest non- 
vanishing) order in tunneling. We calculate the two point 
correlation and four point correlation functions required 
to calculate the cross correlation in terms of edge state 
properties. We analyze the detailed behavior of the cross 
correlation over different Keldysh time domains and pin- 
point the effects of statistics, specifically drawing the con- 
nection between braiding statistics and exclusion statis- 
tics. We show that cross correlation can be expressed 
in terms of two universal scaling functions and that the 
statistical dependence ultimately takes on a remarkably 
simple form. 

It must be remarked that since the perturbation in- 
troduced in Eq. (|3.7|l is relevant in the renormalization 
group sense, perturbation theory is valid only under cer- 
tain conditions. Specifically, perturbation holds as long 
as energies involved, such as voltage and temperature, 
are higher than a cross-over energy scale proportional 



to Fi-^ . For energies much lower than this cross-over 
scale, qp tunneling between edges becomes large enough 
to break the Hall droplet into smaller disconnected re- 
gions and one needs an appropriate dual picture to prop- 
erly address the strong coupling fixed point. For observ- 
ing the effects of statistics as prescribed by the correlation 
calculated here, it is imperative that the system remains 



J 



in the pcrturbativc regime. 

The un-normalizcd cross correlation {Ii){l2)S{t — t') = 
(A/i(i)A/2(t')), where () represents expectation value 
with respect to the full Lagrangian in the presence of 
tunneling events, can be written in terms of expectation 
values with respect to the free theory ()o as 



{h){l2)S{t-t') = -J2{TKhit'^)W"') e'^KL^^*Atl)dt,+^S^L^^,.2{t2)dt2■^^ 



v,v' 



T5](rx/i(i'')e'^-^'"*'^^*^)''*^)o(rA'/2(i'"')e''''-'^-* 



2{t2)dt2\ 



(4.1) 
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By expanding the exponentials in Eq. H4.1|l . we can 
calculate the cross correlation perturbatively in the tun- 
neling amplitude Tj . Since the tunneling operator is rel- 
evant in the RG sense (its scaling dimension is less than 
1), the perturbation theory will break down in the in- 
frared (IR) limit. However, the finite temperature in 
our case provides a natural IR cutoff making this per- 
turvative calculation meaningful. In fact, we will show 



in the following sections that the strongest statistical de- 
pendence is obtained at temperatures comparable to the 
bias voltage. 

Clearly the lowest non- vanishing term in the perturba- 
tion expansion is of order (riF2)^. Using the expression 
for the tunneling current operator in Eq. (|3.8|) . the first 
term of Eq. (|4.1() to this order can be written as 



n-n',vi,V2 «,«' 



(4.2) 
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where we used the fact that the Kcldysh contour time- 
integral in the exponent can be represented using the 
contour branch index rj' 



ri' = ± 



dt 



(4.3) 



and the charge neutrality condition discussed in the Ap- 
pendixiniwhich enforces ei = — e and £2 = — e'- Similarly, 
the expectation value of the current for the second term 
of Eq. (|4.1|) becomes 






(4.4) 



Putting Eq. (|4.2|l and Eq. (|4.4|l together, we arrive at the 
following expression for the lowest non- vanishing term of 
the un-normalized cross current correlation 
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{TKV}^\nV^'Ht'^')vl-'\t'l^)vt''\tT)), 



rie). 



r(-e)/,?;i\ 



Ae'), 



r{-<^') I-lV2\ 



I 



(4.5) 



We now evaluate the two point function and the four 
point function, and then combine them to obtain a sim- 
plified formula for (/i) (12)5(2) (t - f) from Eq. gSJ. 



A. The two point function 

Using Eq. 1)3.6(1 and the result of Appendix [0 to cal- 
culate the chiral boson vertex correlation function, we 
find 



.ieVolXi,*") -ieyo(-^l,tr)\ 



^ {FoF^'r{FoF,-Y'{TK 

X (r^,e"'^i(°^*")e-''^i("^*i')) (4.6) 

= exp[G^,^, (0, t~h)] exp[G^,,,, (0, t~h)] (4.7) 

where G„,„.(x- .t', <- i') = (TKV';(a:,i'')^;(x', f"''))o is 
the Keldysh ordered and regulated ipi propagator for all 
I which will be discussed below in detail for both zero 
temperature and finite temperatures. Here we also used 
the fact that ipi^s are independent from each other for 
different Vs. Notice that the Klein factor contribution 
simply becomes an identity for a two point function of 

I 



the tunneling operator Vi. Hence, Klein factors are not 
necessary in calculations involving two point functions 
alone, and in fact, this is the reason that the lowest or- 
der contribution to shot noise does not contain statistical 
information. However, as we show below, Klein factors 
have non-trivial contributions in the four point function 
which appears in the first term of Eq. ((4.1(1 . 

We now focus on the form of G^_^'(— a,t — t') in the 
limit a ^ 0+. From the detailed analysis of Appendix IbI 
the contour ordered propagators for (jfc and (J)n takes the 
following forms at finite temperatures: 



"sin|[ro + ixr,,^,{t - t'){{t^t')-{x - x')}] 



-i^ln 



{TKC^N{x,f')c^N{x',t'^')) 



/3 



(4.8) 



= lim -i-Xr,,r,'Sgn[wAr(f-t') - (.T-x')]. (4.9) 



Hence the contour ordered propagator for (p represented 
by G in Eq. 1(4.7(1 and Eq. 1(4.18(1 can be written as 



Gr,.r,'{-a,t-t') = -\n 



sinf To + iXr,,r,'{t - t')(t-t' + a)] 



lim «77 ( - + 1 ) X7,.r,'(^ - t') sgn[wjv(^ - t') + 1 



(4.10) 



For the case of interest, i.e., a — > 0+ in the limit tq — ^ 0+, and defined 
the above takes the following form (See Appendix Ib|| : 



G^,^'(0^,t-t') =lnG(t-t'; T,K) 

+ 4x''.v(i-i')sgn(f-f')- (4.11) 

where we have introduced a notation for the amplitude 
of the qp propagator that depends on the temperature T 
and the scaling dimension K = v j]?: 



C{t;T,K) 



(ttto) 



K 



IsinhTrfesTi 



K' 



(4.12) 



X,,,'{t ~t')^[ ^ sgn(t - t') V^ ) . (4.13) 



We can now calculate the expectation value of the tun- 
neling current (/i) of Eq. 1(4.4(1 to order O(r^) using 
Eq. gTI) and Eq. (|^1U|) . 



(/i(f'))(2) =e*|ri|2^e / diie'^""(*-*i)T(0)G(t-ti)2 

(4.14) 



with the statistical factor 



T{9)= [^Tyie'^'^''"' 

\ VI 



,{t-u)sgn(t-t,)\ ^2^sin6' (4.15) 



Notice that the resulting (Ii{f')) is independent of 77 and 
is a constant independent of t, as applicable for a steady 



state current. 



B. The four point function 

Following a procedure similar to the one that led up to 
Eq. 1)4. 7|l . we find that the four point function becomes 



{TkV}'^ it^)%^''^ it"^')vt'^ {tfWt'^ (tf))o 



/y^,g«Vo(^i,t'')g«Vo(^2,t"' )g-iew,(Xi,i;'i)g-«Vo(X2,t^2)\ /2.^g«Vi(0,t'')g-ievi(0,t^i)' 







„jeV2(0,t"' )^-«V2(0,t™)\ 



'l\e' , 



-l\-e'\ 



X {Tkc'' ^^^'^■* ^e-" ^^("''^ ^)o{TKiFoF,-'nFoF,-'Y iFoF,~'r'{FoF^'y ,0 

g[G^,,i (0,t-ti)+G^,^2 (0,t'-t2)+««'G^,^2 (-a,t-t2)+«'G^,^^ (-a,t'-ti)-ee'G^^, (a,t-t')-«'G,i^2 (-a,ti-t2)] 



^2[G^„,(0,t-ti)+G,,,^(0,t'-t2)]. 



.«'[G„„'(-a,t-t')+G'^i„2(-'».*i-*2)] 



(4.16) 



(4.17) 



(TK(Fofri)^(foF2-i)^'(Fo^ri)-^(foF2-i)-^')o (4.18) 



where we used Eq. (IC6|) for the case A^ = 4 to calculate 
the four vertex correlator for the second equality and 
used the fact that Xi — X2 ~ —a. 

The factor involving four sets of Klein factors of the 
form {FqFj'^)'^ in the Eq. H4.18|l has to be appropriately 
rearranged through proper exchange rules, Eqs. (|D2)I and 
(|D5p . as the tunneling operators Vj{t^ys get rearranged 
for Keldysh contour ordering. Since contour ordering 
depends on time arguments of the tunneling operators, 
Klein factors effectively become endowed with dynamics 
in the sense that where each of these four times fall on 
the contour determines whether Klein factors associated 
with tunneling at each time have to be moved or not^. 
In order to rearrange four sets of combined Klein fac- 
tors of the form {FoF~^Y{t^) involved in Eq. H4.18(l . we 
have to take six different possible pairs and contour order 
within each pair. This process can be thought of as a full 
"contraction" where pairs of Klein factors can be treated 
akin to propagators. 

These contour ordered Klein factor "propagators" were 
calculated in Appendix IdI to be 

{TK{FoF-'rif^)iFoF-'r'{t'^'))„ 
{TK{FoF,-'yimFoFr'y'it'^'))o 

_ giee'f sgn(X2-Xi)x^,^/(i-t') _ g-iee'fx,,', „(*'-*) (4.20) 

(4.21) 



where Xv,v'i^) ^^ defined in Eq. (|4.13|) . Further using the 
fact that FqF~^ commutes with itself gives 



{TK{FoFr^y{f'){FoFrY' (t'" ))o = l- (4.22) 

Expressing the Klein factor contribution in Eq. (|4.18l) in 
terms of the six possible pairs and using their forms given 
by Eq. (|4.21|) and Eq. (|4.22|) simplifies it to the form 



-iee'|x„,„2(*-t2)g-iee'fx,i, „'(*!-*') 

-*"'fx^,^'(*-*')e-'"'ix>,i.>,2(*i-*2) ' 



(4.23) 



Hence, we can rewrite the four point function of 
Eq. H4.18|l as the following: 






(4.24) 



e'[G^„,(-a, t-t')-i|x„,^'(t-t')+G„i,„2(-a, ti-ta)-^! X^i ,,,2 (*i -*2)] 



r 



where we used the fact that Gn,ri'{x — x',t — t') = 
Grj',ri{x' — X, t' — t), as it is pointed out in the AppendixlHl 
to replace G,,'_,,j(a, t' — ti) by Gni,r^i{—a, ii — t'). 



Furthermore, observing the following, 



the phase factor by defining the following function 



'2 



e— 

2 



77Sgn(t - 12) - risga{t - t') 
rji sgn(ti - t') - rii{ti - ^2) 



(4.29) 



Grj,r,'{Q , t-t') -■i-^Xv,ri'{t-t') 



-Kin 



sinh ( ii^ 



/3 



tr?j,-?j(0 , t ~ t ), 



+ i-v{l- sgn(i-t')) 



(4.25) 
(4.26) 



V-^i ^&^(^ - *i) + esgn(t - <2) 

- esgn(t - t') - 1} + r/-{l - sgn(t' - ^2)} 

-'7i^{-sgn(i-ii) + esgn(ti - <') 



esgn(ii - ^2) - 1} + m^{'^ + sgn(<' - ^2)}- 



(4.30) 



where we have used Eq. (|4.11(l and Eq. H4.13|l . the four 
point function of Eq. ()4.24|1 further simplifies to the final 
form 




(a) 



(b) 



^ g2[G^^i(0, t-ti)+G^,,^^{0, t'-t2)] 

^ee'lG„,-„iO-, t-t2)+G^i,-„i(0-,ti-t')] 

■^ 

gee'[G,,_„(0-, t-t')+G^i,-,i(0-, ti-t2)] 



sinh 



7r(t-ti) 



/3 



2if 



sinh^^^^ 



2K 



sinh ^^^i^ 



eAT 



sinh^ 



eK 



(4.27) 



f,T7,T7 .TJ1,TJ2 



{t,t',tut2) 



sinh^^ 



e/f 



sinh ^^ „ ^' 



eK' 



(4.28) 



where we defined e = ee' noting that Eq. (|4.24|l depends 
only on the product ee' and collected all contributions to 



FIG. 3: Virtual tunneling events taking place at all times ti 
and i2 allowed by causality contribute to the noise (tunnel- 
ing current cross-correlation) S{t) in the lowest non- vanishing 
order. Solid (red) lines show paths of two quasiparticles tun- 
neling through a FQH liquid from one edge state to two others 
(consitituting /i and I2) and dashed (magenta) lines represent 
virtual tunneling events. Depending on the relative orienta- 
tion of the currents (or equivalently of the virtual processes) 
there are two cases: (a) the case 5* in which two currents are 
in the same orientation (e = -1-1 in the text) and (b) the case 
O in which two currents are in the opposite direction (e = — 1 
in the text). 

Note that e = ee' depends only on the relative direc- 
tions of tunneling, that is e = -1-1 for correlations between 
contributions to /i and I2 in the same orientation (e = e') 
and e ~ —1 for correlations between contributions to /i 
and I2 in the opposite orientation (e = — e'). Each case 
of e = ± is illustrated in Fig. O labeled S/0 respectively. 
The fact that the phase factor $ depends on these rela- 
tive orientations allows us to understand the statistical 



angle dependent contribution to the noise in connection 
with exclusion statistics^ as we will discuss in the next 
section. Another important observation to be made in 
Eq. (|4.3U|I is that the time dependence of the phase fac- 
tor is such that it only depends on the sign of various time 
differences. Hence as we integrate over virtual times ii 
and t2, to evaluate (/i(i)J2(0)) according to Eq. H4.2|l . 
the domain of integration in (ti, t2) space can be divided 
into sub-domains as shown in Fig.^and the phase factors 
for a given set of Keldysh branch indices 77's and relative 
orientation of tunneling e may be computed. Since inte- 
gration is additive, we can consider each domain's con- 
tribution separately and later add them up to unravel 
the effect of the phase factor. A simple analysis to follow 



shows that large parts of the (ti,t2) space do not con- 
tribute to the integral due to vanishing summation over 
Ty's. The remaining domains can then be assigned differ- 
ent values of overall constant from the summation over 
the phase factor. 



C. The normalized noise 

Putting Eq. l|T^ back into Eq. (gSJ for the 0(r?r|) 
term of the current-current correlation and setting t' = 
(the expression has time translational invariancc) and 
taking t > 0, we obtain for (/i(t)/2(0))(2) 



J 






';i.'j2 c,c' 



(t^tq) 



2K 



(t^to) 



2K 



sinhi(^ 



K 



sinh^ 



ifx " 



sinh 



rjt-ti) 



P 



2K 



sinh^ 



2K 



sinh^ 



K 



sinhii^ 



K 



(4.31) 



2e*^|rir2p Yl [dtidt2 e cos[a;o(t - ti - it2)] 



?=±- 



X ( Y^ r]iT]2e 
\'n.';2=± 



j*;''"i'™(t,o,ti,t2) 



(ttto) 



2K 



sinhi(^ 



2K 



sinh^ 



2K 



sinhi^^i^ 



K 



sinh^ 



K\ e 



sinh^ 



K 



sinh 



^{t-t2) 



P 



where we separated the summation over Keldysh branch 
indices y^,„ „i ^ „ from the rest of the integrand us- 
ing the fact that Keldysh branch index dependence en- 
ters the expression only through the total phase fac- 
tor $. To get the second equality, we first used the 
identities ^^ ^,^_^ = J2e e=± (^ = ^^') ^^'^ ^^'^^^ ^^''^ 

Now evaluating the summation over 77's: 



E 

r]i,rj2 = ±l 



mm e 



i$ 



T1,V :i7i,n2 



{tfiMM) 



(4.32) 



for different ranges of (<i, <2),ahows us to extract the 
statistical angle dependence of the cross correlation. 
To do so, we begin by noting that the expression for <1> 
Eq. H4.30|l takes a constant value for different domains 
and hence the sum over 771,772 = ±1 in the brackets of 
Eq. (|4.31|) can be independently evaluated prior to evalu- 
ating the integral if the virtual time (ti, ^2) space is split 
into appropriate sub-domains. Consequently, we make 
the following observations: 



1. We first note that $?' 



^(t, 0,ti,t2) is indepen- 



dent of 772 for <2 > 0. Since I],,2=±i ^2 = 0, 
the expression of Eq. H4.32|) vanishes identically for 
^2 > 0. Hence the only non-zero contribution to 
the integral comes from t2 < 0. 

2. A similar consideration for the summation over rji 
for ti > t allows us to limit the integration range 
to ii < t since ^^''''''''^ •''=' (f , 0, ii X, ts < 0) is inde- 
pendent of 771 . 

3. The above analysis reflects the fact that only vir- 
tual times prior to the times of tunneling events 
can affect the correlation between tunneling events 
due to causality (thus acting as a check that our 
implementation of the Keldysh formalism for Klein 
factors is faithful in respecting causality). 

4. Now we are left only with the region {(ti,t2)| ti < 
t, t2 < 0} which is shaded in the Fig. 0] In this 
region, the phase of Eq. H4.30|l simplifies to 



$;' 



'7,'? ,ni,V2 _/ 



771 1 + ? 



sgn(ii-i2)-sgn(ti) 



-?72 



(4.33) 
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FIG. 4; Subdomains of integration in the virtual time (ti,i2) 
space as it is defined in Eq. I|4.34^ . Only tlie virtual times 
(ii,i2) that lie in the three shaded subdomains Ri, R2, and 
Rs give non-zero contribution to the cross current correlation 
Eq. 14.51 due to causality which is encoded in the Keldysh 
formalism. 



where we used the fact that we are interested in 
the case f > 0. From this expression, which is inde- 
pendent of rj and rj' , we can identify three domains 
that contribute to the integral of Eq. H4.31|l : 



i?l = {(il,t2)Ul<t2<0}, 
i?2 = {(<!, <2)|t2<<l<0}, 
i?3 = {(il,i2)|t2<0<ti<t} 



(4.34) 



as we depicted in Fig. 0] 

(j,^,'?i,j)2 ^^^ Q^ ii, t2) can now be considered as a func- 
tion that depends on the domain Rq {C = 1,2, 3) 
which can be evaluated to give 



$^-'''"-"^[i?i] = %i+r;2), 
$^'''''"-"^[i?3] = 0(r;i+r;2). 



(4.35) 
(4.36) 

(4.37) 



It is worth noting that <I>g has an explicit depen- 
dence upon e only in the domain i?2. However, 
since $g'''^'''^[i?2] becomes independent of 771 for 
e = — 1, the sum Eq. H4.32|l will once again vanish. 
Hence we find that only the case e = +1, where 
the calculated correlation is between the tunneling 
processes with same orientations (case S of Fig.|2Jl, 
contributes to the part of the integral Eq. (|4.31() 
coming from the domain i?2. Later we will discuss 
the implication of this observation in more detail. 

6. Finally we can evaluate the phase factor summa- 
tion over contour branch indices Eq. (|4.32|) for each 
regions using Eqs. (|4. 3 514. 3 7(1 as the following: 



For Ri : J2,j,^,=±iVime 



i^ 



V,V ,^l,T?2 



[-R.1 



For R3 : E^i^,=±i^i'72f 



J2 mV2e''"''^'^'^ = -'isin^ 

?)1,7)2 = ±1 



(4.38) 



?)1,7(2 = ±1 

?)i,r(2 = ±l 



»e(m(i+f)+r,2) ^ J -4sin 6'(2cose') for e = 
I for?=-l 



i6{Vl+n2) 



-4 sin 




r 



At last, we can utilize the observations we made above 
to rewrite the four-point contribution to the noise S{t) 
spelled out in Eq. (|4.31|) in the following form 

(/i(i)/2(0))(^) 
= -8sin2 6'e*^|rir2p V / diidi2? coswo(^-ii-et2) 

-C{h)C{t~t2)'~ 



xiC{t-h)C{t2)r 



C{h^t2)Cit)_ 

16sin^6'cos6'e*^|rir2p / dtidt2cosujo{t - ti - 1-2) 
J R-, 



x{c{t-h)C{t2)Y 



Cih)Cit-t2) 



[C{h-t2)C{t) 



(4.41) 



In Eq. ()4.41f) the statistical angle dependence has been 
extracted out and we found an additional factor of 2 cos 9 
in the contribution from the domain R2 which does not 
occur in the contributions from domains Ri and R3. We 
note that this contribution rises only from the case S of 
Fig. |31 (tunneling processes in the same orientation, i.e. 
e = +1) as it is shown explicitly in the second equality 
of Eq. H4.41|l by the value of e being set to 4-1 in the 
contribution (the second term). 



To complete our analysis, we now turn to the product 
of two point functions (the second term of Eq. (|4.1|) ). 
From Eq. (|4.14() . the second term of Eq. (|4.1() to order 
0(rlTl) becomes 
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(/iW)(/2(o)>(2) 

= -4 sin^ 61 e*^|rir2|^ ^ ee' f dtidt2 
^ ^1 — I- •' 



;ieu)o{t-ti)-ie'u)o(t2) 



e,e'=± 



Cit^tifC{t2f 



(4.42) 



r 



Finally, putting all of the above together, the normal- where SUt) for all values of e = ± and ( = 1,2,3 is 
ized noise (cross current correlation) is found to be of the defined by an integral over domain _R^ for contributions 



following form 

S{i; T/To, K) ^ Y.T. '^c(^~; ^/^o, K) - 1 

+ 2 cos 6* 52^(t; T/Tq,K), (4.43) 



with relative orientation of tunneling e normalized by the 
product of the average tunneling current 



Slii;T/To,K) 



e I diidi2 cos{i - i^ - efa) {C{t - ii)C(t2))' 



C{h)C{t~t2) 



C{ti-t2)C{t) 



E E 



dtidi2 cos[t - ti - l't2]C{t-tifC{i2f 



(4.44) 



e' = ±lC' = l,2,3''-"<' 



r 



where we defined the dimensionless time t = ujQt and also 
introduced a notation for the temperature equivalent to 
the Josephson frequency to be To = fkoQ/kg and used 
the fact that C{t; T, K) = C{i; T/Tq, K) from Eq. ^T^. 
Note that since we chose to study the normalized noise, 
all cutoff dependence and tunneling amplitude depen- 
dence in the noise is eliminated through normalization 
since they enter as common factors for numerator and 
denominator of Eq. (|4.44(l . 



D. Discussion 

Now we can define two scaling functions A and S that 
depend on the scaling dimension K^ the Josephson fre- 
quency Wo = e*V/h (or Tq) and the temperature as 

A{uot; T/To, i^) = E E '^c(^o^; ^/^O' ^) ~ ^ 
S(t; T/To, K) = 2S+{uot- T/Tq, K). (4.45) 

We can thus rewrite Eq. H4.43|l as a sum of two sealing 
functions, corresponding to a "direct term" and an "ex- 
change term" which depends explicitly on cos 6: 



T 
Jo 



cose^fLUot;^,K], (4.46) 



where To = hojo/kB- Hence we were able to extract the 
statistical angle dependence of the noise in the form of 
the factor of cos 9 in the second term while both of the 
scaling functions A and 55 do not depend on the statistics. 
Eq. (|4.46() is the key result of this paper. One can gain 
further insight by recognizing the fact that the domain 
i?i (which contributes to A) and i?2 (which contributes 
to B) are related to each other via exchange of virtual 
times ii and ^2 (See Fig.^. In the following we discuss 
the physics behind the factor cos 6 in Eq. (|4.46|l with an 
eye towards bridging the connection between exclusion 
statistics and exchange statistics in our setup. 

This expression displays a number of noteworthy prop- 
erties: 

a. This finite temperature result is applicable to all 
Jain states (in contrast to recent zero temperature 
results on Laughlin states^SiSS) and it is a universal 
scaling function to the lowest order in perturbation 
theory which is valid provided the tunneling current 
is small compared to the Hall current. 

b. Fractional charge and statistics play fundamentally 
distinct roles, each entering Eq. H4.46|l through dif- 
ferent features: fractional charge through uq = 
e*V/H and fractional statistics through the cos9 
factor. 

c. Given that for Laughlin states 9 < tt/2 and cos 9 > 
0, the exchange term in Eq H4.46|l provides largely 
positive ( "bunching" , boson- like) contributions to 
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the noise whereas for non-Laughhn states with 
9 > 7r/2, its contribution is largely negative ("anti- 
bunching" , fermion-like) . 

d. The fact that only case S, of Fig. [21 contributes 
to this factor can be viewed as a manifestation of 
a generalized exclusion principle^ since the virtual 
processes in this case involves adding a qp to edge 
in the presence of another. It is noteworthy that 
we arrived at an observable consequence of a gen- 
eralized exclusion principle given that Eq. (|4.46() ') 
was derived using anyonic commutation rules pre- 
scribed by braiding statistics^. 

A few remarks are in order here. For the function in 
Eq. (|4.4t)|) to be observable, given that it is the ratio 
between the cross-current correlations and the average 
values of the currents involved, the former needs to be 
at least comparable to the latter. As mentioned above, 
the form of the function only takes into account the low- 
est order terms in tunneling; when higher order terms 
become important, the function no longer remains uni- 
versal in that it exhibits an explicit dependence on tun- 
neling amplitudes. We assumed that tunneling from edge 
occurred from a single point. In our calculations, we 
could have in principle retained the separation scale 'a' 
for points from which qps tunnel into edges 1 and 2, re- 
spectively. Then, 1/a would enter the problem as another 
energy scale. The scaling function would depend on an- 
other parameter ujQa and the cross-current correlations 
would decay with separation length. 

As one might expect from the simple form of the sta- 
tistical angle dependence of Eq. H4.46|l . this expression 
can be understood in terms of simple pictures in which 
the setup simultaneously acts as an accessible stage for 
exchange processes and a testbed for generalized exclu- 
sion. 

To see that the form of Eq. H4.46|l derives directly from 
exchange of qp's between edge 1 and edge 2 via edge 0, we 
first note that information regarding anyonic exchange 
statistics and related time ordering of events is carried 
as a phase factor in the qp propagators. Contour or- 
dering of events and associated time-dependent shuffling 
is completely contained in this exchange statistics phase 
factor. The following analysis of tunneling events and 
their effect on the phase factor not only shows that the 
simple cos 9 factor in Eq. H4.46|l naturally emerges from 
the events associated with exchange statistics, but that 
it is also a manifestation of exclusion statistics. 

Given that causality constrains lowest order virtual 
processes to occur at times ii and ^2 within domains 
Ri, i?2 and R^ of Fig. ^ a qp/qp from edge 1 cannot 
exchange with a qp/qh from edge 2 when ti > (domain 
R3) since qp/qh tunnels between edge 2 and edge via 
virtual process at time ^2 < and back at time before 
anything happens between 1 and 0. Hence it is clear as 
to why S'3 coming from domain R3 contributes to A and 
an uncorrelated two-point piece is subtracted out. 

Now for domains Ri and R2, there is an overlap |ti — ^2! 



in the times during which qp/qh's from edge 1 and edge 
2 stay in edge and hence can experience statistics in- 
duced interference. Between i?i and R2, the only pairs 
of time whose relative orders change from one domain to 
the other is ti and ^2 and hence Ri and R2 are related 
to each other via exchange of entities involved in the vir- 
tual processes occurring at ti and ^2- For both domains, 
depending on the relative orientations of the currents, 
the virtual tunneling processes leave two distinct pairs of 
objects at edge as shown in Fig. 01 

In case S {e = -1-1), where Ii and I2 have the same ori- 
entation , the entities left behind are a pair of two qp's. 
Hence the domain R2 is related to the domain i?i via 
exchange of these two identical particles. This exchange 
between two qp's brings in an additional phase factor of 



in Eq. H4.39(l for the case e 



that adds (con- 



structive interference) to the phase factor e*''i^+"'^^ that 
is common to all three domains (See Eqs. H4.38l4.40|) .'l 
and also to the uncorrelated piece. Since the e indepen- 
dent phase factor sum for domains Ri, —4 sin 9, is the 
same as that of the uncorrelated piece that normalizes the 
noise (and that for the domain -R3), it cancels upon nor- 
malization while the phase factor sum from the case S of 
domain R2 results in the factor of 2cos0 = sin 26*/ sin 6*. 
As a consequence SI is a part of A without statistical an- 
gle dependence while S2 enters 23. In case O (e = — 1), 
on the other hand, tunneling currents have opposite ori- 
entations leaving one qp and one qh behind. Exchange 
between qp and qh brings in an additional phase of e~^^^ 
in i?2 (See Eq. H4.39|l for e = — ) that cancels with the 
common phase factor e"*!^ (destructive interference) re- 
sulting in sin(0 — 6*) = for the phase factor sum. Hence 
there is no 5^ contribution to 15 but only 6*2 constitutes 
S. 

The generalized exclusion statistics is defined by count- 
ing the change in the one particle Hilbert space dimen- 
sions da for particles of species a as particles are added 
(keeping the boundary conditions and the size of the sys- 
tem constant) in the presence of other particlesi^. The 
defining quantity for this description of fractional statis- 
tics is the statistical interaction g^p 



^da^-^g^pl^Np 



(4.47) 



/3 



where {ATV^} is a set of allowed changes of the particle 
numbers at fixed size and boundary conditions boundary 
conditions. For bosons (without hardcore) gap = 0, while 
the Pauli exclusion principle for fermions correspond to 
Qaji = ^a/3- Thus, while the presence of bosons does not 
change the available number of states for an additional 
boson, the presence of a fermion reduces the number of 
available states by one. Once one understands the one 
particle Hilbert space of the particle of interest, the sta- 
tistical interaction gap can be calculated through state 
counting and for FQH qp's such program can be carried 
out using conformal field theory (CFT) as shown by van 
Elburg and Schoutens^. 
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However, the connection between this definition and an 
experimentally measurable quantity has not been clear so 
far. Further, the connection between fractional statistics 
defined through exclusion and anyonic statistics defined 
through exchange has not been clearly established for 
the following reason. The anyonic exchange for FQH 
qp's based on braid groups defined in Refs. [l||3 ] is 
achieved through attaching Chern-Simons flux to hard- 
core bosons^ or spinless fermions^S. However how such 
flux attachments would affect Hilbert space dimensions 
is not obvious. For Laughlin qp's, one can use the knowl- 
edge of the chiral Hilbert space from CFT— to flnd that 
the two deflnitions coincide. Nonetheless, an equivalent 
understanding for non-Lauhglin states are still lacking. 

In our setup, it is quite clear that only virtual processes 
occurring in domains i?i and i?2 with e = +1 ( case S 
in Fig. O host the situation of a qp being added in the 
presence of the other. The difference between i?i and i?2 
is the following: In i?i, a qp from edge 2 tunnels in to 
edge at time i2 and leaves at time with the qp from 
edge 1 present through out. In i?2, however, a qp from 
edge 1 tunnels in to the edge at time ii later "pushing" 
a qp to edge 2 at time <2- Therefore virtual processes 
occurring in domain R2 should be affected by statisti- 
cal interaction of the generalized exclusion principle and 
we can understand why 6*2 alone consists the statistics 
dependent term from this simple reasoning. 

It is astonishingly consistent that the processes con- 
tributing to noise in which the generalized exclusion prin- 
ciple should manifest itself are the ones that are associ- 
ated with the simple factor cos 9 considering that the sta- 
tistical angle is what defines exchange (braiding) statis- 
tics. The fact that our results describe an experimentally 
measurable quantity makes it even more interesting. Al- 
though the manner in which the factor can be related to 
(7a/3 defined in Eq (|4.47() needs further investigation, it is 
clear from our choice of setup which allows access to the 
phase gain upon exchange while satisfying conditions for 
observation of generalized exclusion principle, that the 
connection exists. In the following sections, we analyze 
the scaling functions A{t) and 23 (i) after carrying out 
the integration over respective domains numerically and 
show that the noise experiment of the type proposed here 
can be used to detect and measure fractional statistics. 



THE TIME DEPENDENT NOISE S{t) 
ANALYSIS 



Since the definite double integral in the numerator of 
the expression for Sf Eq. (|4.44() does not grant an an- 
alytic expression in closed form, we evaluate Si numer- 
ically. The integrable singularity in the vicinity of the 
domain boundary is regulated through a short time cut- 
off tq. A natural choice for this cutoff that is intrinsic 
to QH system would be the inverse cyclotron frequency 
uj^^j with ujcfi ^ Imey for the magnetic field of order 
lOT required for the FQH. As we will discuss later, given 



that ojc is such a high energy scale when compared to 
other important energy scales of the setup, namely the 
Josephson frequency ujoft = e*V and the temperature 
ksT, as ought to be the case, all interesting features of 
the noise turn out to be cutoff independent. 

Before we discuss the results of our numerical inte- 
gration, we can gain some insight by looking into the 
asymptotic behavior of SUujQt) in the short time limit 
ksTt, ojot <C 1 and in the long time limit ksTt, uJot » 1. 
In the short time limit, the numerator of Eq. 14.44|l be- 
comes 



Sl{t)(xi\kBTt\'^ / dtidt2Cos{ti+it2 



where we used 



Re 



lim C(t) ex InkBTt 



C{h)C{t2) 



Citi-t2) 



-~K 



(5.1) 



(5.2) 



Since the domain size of R^ is proportional to t, the inte- 
gral in Eq. \b.l\ is negligible for C = 3 and hence we can 
ignore 5*1 's contribution to yi as t ^ 0. However since 
the integral in Eq. (|5.1|) is independent of t for domains 
Ri and i?2, Si{t) diverges as (ksTt)^^ for e = — as 
i ^ and both S^{t) and 5^ vanish as (ksTt)^ . As a 
result, A{t) defined in Eq. H4.45|l is dominated by S^ in 
the short time limit and we find 



(5.3) 



\A{t)\ ^\kBTt\~'' il + ---). 

On the other hand, 

mt)] = 2\S+it)\ ^ IfcsTil^. (5.4) 

In the long time limit, the numerator of Eq. H4.44|l can 
be approximated for C = 1, 2 {SI and 5^) by 



SUt) 



~ „~2Kir\kBTt\ 



diidt2 cos(< — fi — £^2) 



Ri 



C{i2? 



C{h) 



C{tl-t2) 



(5.5) 



which vanishes as e-'^K'^\kBTt\ niultiplied by an oscillat- 
ing factor. Also, since the size of the domain R^ grows 
with f, Sg and S^ become finite and we find 



yiu,[s+{t) + s^[t)]-i = Q. 

t — 'OO 



(5.6) 



Therefore one can expect both A{t) and 23 (t) to vanish 
exponentially in the long time limit due to thermal fluc- 
tuations. 

Depending on the relative magnitude of the two energy 
scales of the problem, namely the Josephson frequency 
Lo^h ~ e*V and temperature kgT, there may or may not 
be an intermediate range of times over which both A(t) 
and 23(t) show oscillatory behavior. Figure O shows the 
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A{iJot), -Biujot) 



0.25f 



-0.25 



-0.5 




(a) High temperatures. 



A{ujot), 'B{uJot) 






Loot 






(b) Low temperatures. 



FIG. 5: The scaling functions A(ujot) (blue dashed line) and 
'B(ajoi) (green solid line) for v = 1/5 as a function of dimen- 
sionless time cjot at two different temperatures: (a) T/To = 2 
and (b)T/ro = 0.25 where the reference scale is set by the 
Josephson frequency to be To = ^^. A{iJot) is dominantly 
negative at short times due to the divergent contribution com- 
ing from counter-current correlations while ^{ujot) coming 
only from correlations of currents in the same direction is 
dominantly positive at short times. A{ujat) and ']i{ujot) for 
other filling factors show similar forms up to an overall scale 
set by the associated scaling dimension K. 



'B(t) are of comparable magnitude hence the cos 9 fac- 
tor in the "exchange term" of Eq. H4.46|l can significantly 
shift the total noise positively or negatively depending 
on whether 6 < tt/2 or 6 > tt/2. Hence it is important 
to understand the intermediate time scale and it can be 
done so numerically. 

For temperatures T <^ Tq — ujQh/kB, both A{t) and 
'B(i) in Fig. [S] (b) demonstrate definite oscillatory be- 
havior characterized by the Josephson frequency ujq with 
an exponentially decaying envelope q-K'^^bTi ^ Although 
the magnitude of !B(t) is much smaller than that of yi(t) 
due to the dominant short-time divergent contributions 
to A{t) from S^ , one can expect the Fourier spectrum of 
A{lo) and 23(0;) to show singular behavior in the vicinity 
of cj = wq. Such features of the Fourier spectrum can be 
used in experiments to locate the Josephson frequency 
which is determined by the fractional charge of the FQH 
qp. In the next section, we discuss the Fourier spectra 
A{uj) and "^{oj) with a particular interest in extracting 
signatures of fractional statistics and fractional charge. 



VI. FREQUENCY SPECTRUM S{lo) 

As with the time dependent noise, it can be seen from 
Eq. H4.43|l that the frequency spectrum of the normalized 
noise also gets contributions SI from different domains 
i?^ as the following: 

/oo 
dte"^^S{t)= Y^ Ysi{u) + 2coseS+{Lu). 

(6.1) 
Noting the fact that the oscillatory time dependence of 
all Sl{t) comes from the factor cos[wo(i — ii — 6^2)] in 
Eq. H4.44|l [t = ujQt), we can use the identity 



typical time dependence oi A{ujot) and 'B(woi) as a func- 
tion of dimensionless time ujQt = t ai high temperatures 
LUoH ~ ksT and low temperatures ujoH 2> ksT. 

We can confirm from Fig. \E\ (a) that when the charac- 
teristic time for the exponential decay in the long time 
limit {2KTTkBT)~^ (See Eq. (|5.5|l ) is of the same order or 
less than the oscillation period o;,^^, no oscillations are 
seen. At the same time, the short time behavior which 
sets in only for t ^ {ksT)^^ is a less prominent feature 
of the plot. At moderate time scales however, A{t) and 



cos[a;o(t — ti — £^2)] = cos(woi) cos[wo(^i + ^^2)] 

+ sin(a;oi) sin[wo(ii + £^2)] (6.2) 

to factor out oscillatory factors from each SUt) as the 
following: 

Sl{t) oc [cos(wot)5'c(i) + sin(cjo<)3<c(0] , (6-3) 

where the decayin envelope functions 3^5 {t) and ?{J (t) are 
defined by 



Jlit) = I / dtidt2 cos[wo(ti + et2)] {C{i- h)C{t2)y 
:Kl{t) = e [ dtidt2 sin[u;o(ii + £^2)] {C{i - h)C{i2)f 

jRc 



C{h)C{t~i2) 



C{ti^t2)C{t) 
C{h)C{i-i2)' 



C{h~t2)C{t) 



(6.4) 
(6.5) 
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from Eq. |(n3 and Eq. (|0^ . 

It is clear from Eq. (|6.3(l that 5(w)5 is a convolution 
of an oscillating part and the envelopes which we may 
simply compute. Thus we have 



5H^cx- r^{u + loo) + 9liLu - too) 



2i 



:kI{lo + loo) + :KliLo ^ loo) , (6.6) 



with 

3<l{uj) 



dtTAt)i 



= 2 



dt cos Lot 'J% [t) 



dt D<1 (t) e'"* ^2i dt sin Lot J{J (t) , 



(6.7) 



where we used the fact that 5"^ (t) is even in t while IK^ (t) 
is odd in t. Hence we can understand the frequency spec- 
trum by investigating frequency spectra 3^Ulo), !K5(w) 

I 



and superposing them with their centers shifted by itwo- 
One important observation to be made in Eq. (|6.7|l is the 

fact that S^Uoj) is an even function of w while 'KUlo) is an 
odd function of ti;. As we discussed in the previous section 
(See Eq. (|531)), S'^(i) decays exponentially as g-^^'^^'^s^* 
in the long time limit {KnkBTt) ^ 1 and hence 



Flit) ~ :Km 



o 



'2KnkBTt 



(6.8) 



in the long time limit as well. As a result, the form of 
3^1 (i^) is a broadened peak at w = with the width of the 
peak 2K'KkBT determined by the temperature and the 
scaling dimension of FQH qp. Similarly, !K^(w) takes the 
form of a derivative of a peak with the same width. The 
characteristic behavior of 9^5 (a;) and "K^Ilo) is plotted in 
the Figure [3 

From the definition of the "direct" and "exchange" 
terms in Eq. (|4.45|) . we can express A{lo) and 23 (w) in 
terms of 5'J(w) and "KUlo) as 



C=l,3?=± 



s 



^=1,3 e=± 

(o.) = i [j+(a.o + c^; T, A') + ^+{uoo - ^; T, K) 

^+{uoo+uj- T,K) - S+(^o - uo- T,K) 
I 



1 
2i 



(6.9) 



(6.10) 



In the rest of this section, we use the properties of !?^(i^) 

and "KUijo) to understand and discuss A{lo) and 'B(w) in 
two frequency regimes of interest: low frequency lo -^ loo 
and near Josephson frequency w ~ wq- 

For w ^ Wo, both yi(w) and Ti{ijo) collect tail ends of 
shifted 3^J(w) and "KUlo). Naturally, they show white 
noise behavior and hence we focus on zero frequency val- 
ues yi(a; = 0) = Ao and ®(a; = 0) = !Bo which depend only 
on the temperature(See Fig. 0). The direct term ^Iq is 
negative due to its dominant contributions coming from 
e = — processes in domain Ri with the opposite orien- 
tation of the currents as illustrated in case O of Fig. O 
On the other hand, 23o is positive as it only involves case 
5(e = +) from domain i?2. As a result, Eq. (|4.46|l trans- 
lates into 



5(^ = 0) = -|yio|+ cos 6*1201 



(6.11) 



for the frequency spectrum of the noise at zero frequency. 
At low temperatures (T < To, To ee fiwo/^s ^ SOmif 
for F = 40/iy), |yio| dominates the spectrum since con- 



tributions from HCf (wo) and ?']~(wo) both diverges in the 
T ^ limit displaying T~^ power law divergence. How- 
ever at higher temperatures, both 3'^(w) and !K5(w) are 
more like Lorcnzian for all e and C,. Hence for tempera- 
tures of order Tq, So is comparable to yio as we display 
in Fig. |S1 (At even higher temperatures, both correla- 
tions are washed out by thermal fluctuations.) Conse- 
quently the statistics dependent exchange contribution, 
cos 6* I So I, will noticeably affect the experimental mea- 
surement of the total noise at small frequencies at such 
temperatures. It is rather astonishing that the statisti- 
cal angle 0, which was originally quantum mechanically 
defined in a highly theoretical basis, can be measured in 
this reasonably accessible frequency range at finite tem- 
perature. 

Near the Josephson frequency, A and CB develop qual- 
itatively different sharp features at temperatures T <^ 
To/K: A{uj) nearly crosses uj axis steeply while "BIlo) 
develops a peak of width 2KkBT/h. This can be un- 
derstood by analyzing which terms among 5"^ and 3^5 
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^0, 



T/To 



FIG. 6: A typical temperature dependence of zero frequency 
values of the direct term 71 (lj = 0) (orange dashed line) and the 
exchange term B (c^j = 0) (green solid line) taken from v = 2/5 
case. The exchange contribution 2 (tj = 0) is comparable to 
the direct term in the temperature range 2 < T /To < 4 and 
hence the fractional statistics has a marked effect on the total 
noise at small frequencies in the same temperature range. The 
qualitative behavior is universal to all fractions. 



dominate each oiA{uj) and ^(a;). As for A^ its dominant 
contribution is from domain Ri with e = — . Further, 
i!Kj~(a;) is strongly dominant over J^{uj) in magnitude 
as can be reasoned as below. _ _ 

Recalling from Eq. (|6.5|l that 3'J~(w) and iJC^(a;) 
have cos ((jJo(^i ^ ^2)) \svah{TTkBT{ti — ^2))!" 

and sin (cjo(^i — ^2)) |sinh(7rfcBT(ti — ^2))! 

factors in their integrands respectively, 
of their characteristic behaviors can be 

h 



tured 



by 



dt'cost'|sinh(7rfcsri') 



each 

cap- 

and 

domi- 



J„ dt' smt' \smh{TikBTt')\ with the latter 

nating over the former for {KTrksT) <C 1. 

Comparison between 3^^(w) and i'yC[{uj) plotted in 
Fig. [7| this observation and hence A{uj) in the vicinity 
of oj = tJo can be approximated as 



A{Lo)K-'K^{uja 



CO) 



(6.12) 



plus small corrections from other contributions including 
the tail of the same function centered around lu ~ —ujq: 
5{j~(wo + 1^)- Since OiUoj) is an odd function of lu as 
we mentioned earlier, the dominant contribution from 
!Ki {ujQ — Lo) to A{lu) makes it nearly cross uj axis in the 
vicinity of the Josephson frequency and thus suppressing 
the contribution oi A{u}o) to the total noise. 

We now turn to 'B which determines the mag- 
nitude of the "exchange" term that bears the sta- 
tistical angle dependence. Given that only do- 
main i?2, e = + contributes to 23, the parts of 
the integrands that strictly depend only on ti and 
t2 now becomes cos (^0(^1 + ^2)) |sinh(7rfcBT(ii — ^2))! 
^2)) |sinh(7rfc_BT(fi — <2))| respectively 
Since the argument of the trigono- 
metric function is the sum {ti + t^) while the difference 
(ti — ^2) enters sinh function in a non-singular fashion, 
both 3^2 (<^) 3'iid !K2 (^) arc of comparable magnitude 
away from their centers uj = (where 3^2 (^ = 0) 7^ 



and sin (cjo(ii + 
for ?+ and 3^+ 
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FIG. 7: Two dominant contributions to yi as a function of 
uj/uiQ at temperature T = O.STq: 5'^(aj) in (green) solid line 
and i'K^^{uj/uJo) in (blue) dotted line. The width of the peaks 
are 2KkBT. i^^ (uj) is already dominant with stronger con- 
tribution at the tail even at this temperature. As the tem- 
perature is further lowered, the dominance of iJi^ (uj) over 
3^^ (io) becomes even more prominent. 



(even function) while !Hj(w = 0) = (odd function)). 
Thus 23 (w ~ Wo) can be appropriately approximated as 



'B{uj) «3^+(cjo-cj) 



(6.13) 



which has a peak centered at lu = luq- 

From the preceding analysis, we can conclude that the 
exchange term plays a significant role in the measured 
net correlation S(lu/luq = 1) by shifting it to a positive 
(negative) value for Laughlin (non-Laughlin) states since 
the exchange term has cos 6' > (cos 6' < 0) and the 
direct term A{lu) is strongly suppressed near Josephson 
frequency. These sharp features not only provide a di- 
rect measurement of the fractional charge associated with 
the Josephson frequency but also allow a natural con- 
trast between two filling factors with the same charge 
but different statistics, e.g. v — 1/5 and u — 2/5 
as we show in Fig. |HIa). The information on statistics 
and charge contained in the noise near Josephson fre- 
quency is recapitulated in Fig. |HIb) by introducing an 
effective charge e = fvJj/V from the crossing point of S 
i.e. S{lu/luq) = 0). The deviation between e for i/ = l/5 
and ly = 2/5 satisfying ei/5 < e* < 62/5 over a broad tem- 
perature range despite their common fractional charge 
e* = e/5 is a consequence of different statistical angles: 
9 = 7r/5 (bunching) for 1/ = 1/5 and 6 = 3tt/5 (anti- 
bunching) for v~2/b. Observation of such deviation will 
be a definite signature for the existence of anyons. 

We conclude this section by summing up the two main 
results which show ways to extract the effect of fractional 
statistics: S{lu = Q) = — l^lol + cos6'|'Bo| for the zero fre- 
quency Fourier spectrum plotted in Fig. [^ e* defined 
from S{ijj) for lu near Josephson frequency. The first re- 
sult can be used to measure the statistical angle in a 
quantitative manner while the second result can be used 
to determine the sgn(cos 6) and hence detect bunching 
versus anti-bunching behavior. 
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FIG. 8: (a) S{uj/loo) for u = 1/5 and ;/ = 2/5 for T = 
20107^,40017^ respectively with V ~ 40/xF. S{ijj/ujo) shows 
a peak in the vicinity of the Josephson frequency uoo ■ At u)o , 
S{ui/ujo = 1) is positive for the Laughlin state v = 1/5 and 
negative for the Jain state u = 2/5. (b) The effective charge 
e(r) in units of e is determined by the crossing points of (a) 
for 1^ = 1/5 and i^ = 2/5. 



VII. CONCLUDING REMARKS AND 
MULTI-TERMINAL NOISE EXPERIMENTS 



In summary, we have shown that the cross-current 
noise in a T-junction can be used to directly measure 
the statistical angle 9 independently from the fractional 
charge, and conceptually, from the filling factor v. The 
normalized cross-current noise as a function of frequency 
oj assumes the form of a scaling function that depends 
on the quantities uj/loq and HiOo/kT, where 'wq' is the 
Josephson frequency, and all the statistical dependence 
is contained in a simple cos 9 factor. The low fre- 
quency results show that there is a quantitatively sig- 
nificant statistics dependence of this noise. We have also 
shown how to measure both fractional charge and statis- 
tics near the Josephson frequency. Specifically, since 
the Josephson frequency explicitly depends on the frac- 
tional charge, identifying its approximate value based 
on where the finite frequency cross-current noise drops 
to zero would provide a means of measuring fractional 
charge. It must be noted that the fractional charge, being 
a single-particle property, can be detected through auto- 
current correlations, i.e., shot noise measurements. The 
zero frequency shot noise is proportional to the fractional 



charge in the weak tunneling limilsiSiiSiii^ while the fi- 
nite frequency shot noise shows a peak at the Josephson 
frequencjiS^. 

Plotting the actual frequency at which the cross- 
current noise drops to zero as a function of temperature 
and comparing this value to the expected Josephson fre- 
quency would give a measure of the bunching versus anti- 
bunching behavior of the qp's arising from statistics. A 
definitive means of measuring such statistical behavior 
would be by comparing the cross-current noise of two 
FQH states with the same fractional charge but different 
fractional statistics, thus pinpointing the effect of topo- 
logical order. We remark that our results are independent 
of non-universal short time physics. 

In this paper, we have considered a setup, the T- 
junction, consisting of quantum Hall bar with three 
terminals. In this setup, two quasiparticles (or a 
quasiparticle-quasihole pair) are created (or destroyed) 
at nearby places of one terminal (terminal in Fig. 2). 
In the final state the two other terminals (1 and 2 in 
Fig. 2), which are widely separated, each detects just 
one quasiparticle. We have assumed throughout that the 
two quasiparticles are created at nearby points in ter- 
minal separated by a distance a small compared with 
the magnetic length. Our calculation can be straightfor- 
wardly generalized for the situation in which a is some- 
what greater than the magnetic length. However, as 
a increases, various decoherence effects will hinder the 
detection of the statistical phase. In particular, when 
a > vH/kBT where v is the speed of the propagating 
charge mode, thermal fluctuation will overwhelm the co- 
herent propagation of edge state between two points, and 
the interference effects associated with fractional statis- 
tics will become rapidly unobservable. 

Safi et.al^ have recently considered the identical setup 
for Laughlin states at T = and obtained a series expres- 
sion for S{lli ~ 0) for the case O alone as a function of i/ 
(note that for Laughlin states one cannot distinguish v 
from K or e*/e or O/n). Indeed, our finite temperature 
calculation shows that the total cross current noise will be 
completely dominate d by case O as T ^ which was the 
case of interest in Ref. |23 . However, our study also shows 
that, at finite temperatures, case S brings in an explicit 
statistics dependence in the form of cos 9 which can be 
distinguished from the effect of the fractional charge and 
the scaling dimension for non-Laughlin states. Further- 
more, even at relatively low temperature at which most 
noise experiments involving QH systems take place (see 
below), this statistics dependent contribution was found 
to be of comparable magnitude as the statistics indepen- 
dent contribution. It is rather encouraging to find that 
a purely quantum mechanical property such as statistics 
can in principle be observed in our setup at non-zero 
temperatures. In fact it is precisely the finite tempera- 
ture that assists processes that are sensitive to statistics. 
A closely related setup is the four terminal case studied 
(onl y fo r Laughlin states at T = 0) by Vishveshwara in 
Ref. [23 where the cross current noise at equal times t — 
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was calculated for the contribution that involves all four 
terminals. We should note that however, even in this 
case, an actual measurement of the cross-correlation will 
have contributions from correlation between three of the 
four terminals which is precisely S(t) we have calculated 
here. A number of interesting interferometers have been 
proposed to detect fractional statistic32Si2i2^i2^, none of 
which has yet been realized experiment ally, ex cept pos- 
sibly for the recent experiments of Ref . [3al37ll3a | . The- 
oretical schemes to measure non-Abelian statistics have 
also been proposed i^2ii£ 

Here we have presented the calculation using the edge 
state theory constructed by Lopez and Fradkin for Jain 
stateaii. However, as long as one can limit oneself to a 
single propagating mode by either choosing the most rele- 
vant mode of hierarchical picture (see below) or by using 
Lopez-Fradkin picture which only has one propagating 
mode by design, the theory presented here is applicable 
and could indeed be used to compare predictions from 
different descriptions. 

For the T-junction setup studied here, the hierarchical 
picture would imply multiple prop agating edge modes 
for non-Laughlin states. In Ref.[l3, Wen derived the ef- 
fective theory for edge states of hierarchical FQH states 
following Haldanes hierarchical construction^. Accord- 
ing to the hierarchical theory, the i/ = 2/5 FQH state is 
generated by the condensation of quasiparticles on top 
of the i^ = 1/3 FQH state. Thus the 2/5 state con- 
tains two components of incompressible fluids. If we 
consider a special edge potential such that the FQH 
state consists of two droplets, then in this picture, one is 
the electron condensate with a filling fraction i^i = 1/3 
and radius ri, and the other is the quasiparticle con- 
densate (on top of the 1/3 state) with a filling frac- 
tion t^2 = 1/15 ( note 1/3 -|- 1/15 = 2/5) and radius 
r2 < ri. Recently, one of us used this picture^ to ac- 
count for the superperiod Aharonov-Bohm oscillations 
observed in recent experiment32SiSi2& where the smooth 
edge potential could indeed have created two droplet sit- 
uation. On the other hand, if the edge is sharply defined 
so that there can be direct tunneling between descendant 
states (2/5 state for example), it is reasonable to consider 
the tunneling current being dominantly carried by the 
mode which has the smallest scaling dimension (that is 
"most" relevant in the renormalization group sense). In 
Wen's theory, although the primary edge qp mode carries 
e* = 1/5, 6'/7r = A' = 3/5, it is actually the qp mode with 
e* = 9/-K = K = 2/5 that is most relevanliiii, which is quite 
different from the e* = 1/5, 6'/7r = 3/5, AT = 1/10 for the 
charge mode in Lopez-Fradkin picture. In fact, given that 
our calculation is applicable to alternate pictures which 
come with different predictions, an experiment such as 
the one proposed here, where the charge and statistics 
of Jain state quasiparticles can seperately be identified, 
becomes all the more immediate. 

Turning to experiments, some recent mesurements 
in a quantum interferometer geometry were interpreted 
to imply fractional statistics. An effective charge ex- 



tracted in recent shot noise experiments— was somewhat 
larger than the minimum quasiparticle charge, suggesting 
bunching behavior. Observation of Aharonov-Bohm os- 
cillations with a superperiod in an anti-dot ^Si^L2& have 
also been ascribed to statistical effects. However, the 
noise measurements we propose here will allow for a di- 
rect and independent way of measuring both fractional 
statistics and charge separately. Low frequency predic- 
tion have the advantage of yielding quantitatively signifi- 
cant statistics dependence of the noise. Near the Joseph- 
son frequency, our prediction has the strength of probing 
both charge and statistics simultaneously. The geome- 
try proposed and the required experimental conditions 
are all very much within current capabilities. The T- 
junction setup is reminiscent of the multi-lead geometry 
employed in effective charge measurement, though, the 
issue of whether or not the fabricated T-junction would 
form a single point contact would be relevant. Typical 
parameters used in noise experimentsi^i^ are bias volt- 
ages of 10 — 150/iy (corresponding to a Josephson fre- 
quency ojo ~ 1 GHz), and temperatures of 10 — lOOmK. 
These parameters access the ratios of T/Tq shown here, 
making our proposal plausible. The parameter regime in 
which the perturbative results for cross-current correla- 
tions are valid depends on the inter-edge tunneling am- 
plitude and consequently on fabrication. The condition 
(on applied voltage and temperature) for remaining in 
this regime is that the ratio between tunneling current 
and maximum edge state current ve^Vjh be small, i.e., 
of order 0.1. 

A definitive measurement of anyonic statistics, such as 
the one proposed here, would be one of the first instances 
of experimentally establishing the existence of topologi- 
cal orderiiSiiS. Employing topological order has recently 
emerged as a possible route to decoherence free quan- 
tum computing^. Clear signatures of fractional statistics 
would mark the very first steps towards bringing some of 
these fascinating ideas to life. 
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APPENDIX A: CHIRAL AND NON-CHIRAL 

BOSON COMMUTATION RELATIONS AND 

CORRELATION FUNCTIONS 

In this appendix we quantize the chiral boson theory 
and calculate its correlation function. We also extend 
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this calculation to the case of a neutral topological field 
that is non-propagating. 

The Lagrangian density for a left/right moving free 
massless chiral boson field (f)± (x ± vt) propagating with 
speed V is 

C± = -^dA±{±dtd^± ~ vd^^±) (Al) 

47r 

with the corresponding action Sq — J cPrC where f = 
(x, t) and g is a positive real parameter. Notice that the 
Lagrangian is invariant under an arbitrary overall shift 
of the field (f>±, i.e. 



>C±[</'±(r)]=/:±[(/)±(f) + «] 



(A2) 



with an arbitrary constant a. Hence C± is independent 
of the zero mode of cl)± . 

Correct quantization of these chiral bosons can be 
achieved by the following equal time commutation re- 
lations 

[0±(x,t),±|-a;</>±(x',t)] = liS{x~x') (A3) 

where the factor of ^ comes from the fact that (j)± 
are non-local fields with respect to the non-chiral bo- 
son = 0_ -I- (/)_!_ occupying only half the phase space 
of (f). The non-chiral boson is a true local field and it 
can be quantized in the usual way. We will find below 
that Eq. (|A3|1 is indeed the correct commutation relations 
consistent with the commutation relations of non-chiral 
boson. To find the non-chiral boson Lagrangian using a 
path integral, we need another (dual) non-chiral bosonic 
field 9 to change the basis from </)_ and 0-|_ to (f> and 6. 
Hence we define 



(A4) 



and write (j)± in terms of (f) and 6 as the following. 



0+ = i(0-0), 0_ = i(<^ + 0) (A5) 

Then the total Lagrangian for two fields can be written 
in terms of and 9 as the following 

/:+ + £- = -^[d^(t)dt9 + d^9dt(i) + v{d^<i)f + v{d^ef] 

.(A6) 
and the path-integral for two chiral fields can be written 
as a path-integral for two non-chiral fields to give 



D0_ / T»0+e'i''''''[^++^-l 



D0 / D6)e~'itJ'*''-[''(9^e)'+29*</'ax«+i'(ax0)"]^ 



(A7) 
Now we can integrate out from the above equation 



to obtain the Lagrangian for the non-chiral field 



(A9) 



We can find the canonical conjugate of 4>, H^ in the usual 
way to be 



H. ^ 4^ = if a.0 



(AlO) 



with the equal time commutation relation 

[cbix),n^{x')]=i5{x-x'). (All) 

However, since (j) = (j)-{x — vt) + 0+(a; + vt) 

dt<p = dt{<p- + <p+) 

= -vd,{(f)- + cf)+) (A12) 

= —vdx9 

which means dx9 is proportional to the canonical conju- 
gate H^: 



H, 



-i"' 



(A13) 



and hence Eq. IjAlip can be rewritten as 



47r 
[(pix), dx'9{x')] == -~i—S{x - x'). (A14) 

g 

From Eq. IjASp and Eq. (|A14|) one can check that Eq. (|A3p 
indeed is the correct way of quantizing the chiral modes. 
Now integrating both sides of Eq. HA3|I with respect 
to x' — x and noting that dx'-x[4'±ix)T4'±i^')] = 
[(j)±{x) , dx' <l)±{x')] — [dx4>±{x) , (l)±{x')], one finds the fol- 
lowing equal time commutation relations for left/right 
moving bosons 4>± 



[(j>+{x),(l)+{x')] = -i-sgnix-x') 
9 

7T 

[0_ (x), 4>_ {x')] ^ i- sgn(x - x'). 
.9 



(A15) 



To calculate the chiral boson propagator, we first de- 
fine a generating functional: 



Z[a] = /X'(/)± e 



J J d^r[C±[ct>±]+a(j}± 



Z[0] e 



J d'^r f d^r'a{r)M^^(,r-r')a{r') 



(A16) 



where we defined the differential operator associated with 
the Lagrangian C± 

9 



M± = ij-dx{±dt-vdx). 



(A17) 
Now the propagator is nothing but the inverse of M since 



G±(f-7^) = (0±(r')0±(r')) 
5 In Z\a\ 



5a{r)5a{f') 
Mi\T-^), 



(A18) 



a=0 
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and can be calculated using a Fourier transformation 
with the proper epsilon prescription for a time ordered 
propagator (substituting oj with oj + ie sgn(w) in the limit 
e->0+): 



G±{r) 



dkduj 2n 



^i{kx—ujt) 



(27r)2 ig ik[^i{uj + iesgn(a;)) — ivk] 

dkduj 1 e-tiu;t±kx) 



lix g k(uj + ie sgn(tj) — vk) 
k 



(A19) 



Here we introduced a UV (short distance) cutoff tq and 



used the fact that the lu contour integral has to be closed 
in the lower(upper) half plane for i > (i < 0) and 
sgn(a;) = sgn(fc) near the pole at w = vk. Notice that 
the left (right) moving propagator G+ {G-) indeed de- 
pends only on vf+x [vt—x)^ i.e. 4>±{x,t) ~ (j)±{vt±x). It 
is clear from Eq. HA19|I that G{r) has a logarithmic diver- 
gence in the small k (IR) limit. However, one can define 
a finite quantity by subtracting the equal position propa- 
gator G(0), which also has a logarithmic divergence, and 
obtain the following 



G±(r) = G±(r)-G±(0) 
1 f dk 



9 



1 



■In 



[{e,-^k(vt±x)-kn-, _ i}e(fc)e(^) _ ^^^k(vt±x)+kro _ i}e(_fc)e(_t)] 

To + isgn{i){vt ± x) 
To 



(A20) 



r 



At finite temperatures T, the regulated and time or- 
dered propagator can be calculated from Eq. ljA20p via 
conformal mapping to take the following form 



G±{x',t)^--\n 
9 



■pTQ 



sm 



llf+^sgn{t^{vt±x)^ 

(A21) 
where (3 ~ l/{kBT) is the inverse temperature with ks 
being the Boltzmann constant. 

In the non-propagating limit of vm -^ 0+, Eq. (|Aip 
becomes the Lagrangian for a topological fieldiL^ as the 
following 



the chiral field: 

i r dkduj e^(kx-ujt] 
gN J 27r kuj 

= — - — {~iTTQ{t) + iTTQ{-t)){i-K<d{x) - iTTQ{-x)) 
9n Stt 

% 7T 

= 77 sgn(i) sgn(a;) 

gN ^ 



lim — -sgn(t)sgn(u7vt + x), 
'«jv^o+ gj^ I 



(A23) 



where the limit should be taken at the end of any calcu- 
lation that uses ((/)Ar(r)0jv(O)). 



APPENDIX B: CONTOUR ORDERED 

(KELDYSH) PROPAGATORS FOR CHIRAL 

BOSONS 



r 9n ^ 

l-N — -^Ox(, 



>N(Jt(pN 



(A22) 



where g^ can be either positive or negative real number. 
Cf4 is still invariant under an arbitrary overall shift of 
the topological field (/)jv and the propagator for the topo- 
logical field can be calculated in the similar way as for 



Here we derive the contour ordered Green functions 
for chiral bosons described by the Lagrangian Eq. (|A1() 
and topological modes described by Eq. HA22|) . Time 
ordered propagator can be thought of as a sum of two 
terms depending on the sign oi (t — t'): 

{Tcp{x,t)<p{x' ,t')) 

= Q{t^t'){^{x, t)${x', t')) + eit'-t){${x', t')4>{x, t)). 

(Bl) 
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In the same way, a contour ordered propagator is deter- 
mined by the order in which different times appear on 
the contour sketched in Fig. El 



forward rj — + 

>- 



-«- 



J 



backward r/— 



FIG. 9: The Keldysh contour from t = 
forward branch (?? = +), and back to t 
branch (?; = —). 



-co to i = CO along 
-co along backward 



G^,n'{x^x',t-t') = (rK0(x,f' )0(^',i"'')) 



{<l)ix,f')q^{x',t"'')) fori" >Ki"'' 
{(j){x',t"^')(bix,f')) fort'"' >A'i" 



(B2) 



where we used the expression t >k t' to denote t appear- 
ing later than t' in the contour. Clearly when both times 
occur in the forward branch (77 = r/' = +). contour or- 
dering is the same as time ordering and when both times 
occur in the backward branch [r] = rj' = —) , contour or- 
dering is the same as anti-time ordering. On the other 
hand, when t is on the backward branch while t' is on the 
forward branch {tj = —, rj' = +), obviously t is occurring 
later in the contour and the exact opposite is the case for 



r] = +, rj' = —. Comparing the time ordered propagator 
Eq. I|A20(I with Eq. ||BT|| and setting v = 1 



l±{x,t)^±{0,0)) = --\n 
9 



To + i{t± x) 



To 



(B3) 



and combining this with the above observation, we arrive 
at the following contour ordered propagators for the right 
moving chiral boson (/)„ : 



J 



G-+(x,t) ^ (0_(a;,f)0_(O,O)) = - In 

.9 



G+-(a;,f) = 


= (0_(O,O)(/)_(a;,<)) = -In 
9 


G++(a;,f) = 


= (T0_(a;,f)0_(O,O)) = il 
9 


G..{x,t)-- 


= (T</)_(0,0)^_(a;,f))-il 



To + i(t- 


-X) 


To 




To - i[t - 


-x) 



In 



In 



To 



To -f isgn(t - 0)(f - x) 



To 



To 



(B4) 



To - isgn(t)(f -x) 



which can be written in the following compact form 



G^,,,'(x,t) = — In 
where we introduced 



To +ix»;,r,'(0(^-2:) 



To 



(B5) 



Xri.ri'it) = -——Sgn{t) -— 



Gn,n' {x, t) = lim -i-Xv,v' (*) sgn(uArt - x) (B7) 



(B6) Notice that both G_+ (a;, 0)-G+-(a;,0) and G^+(a;, 0)- 
G^_{x,0) properly yield appropriate commutation rela- 
Going through the similar analysis for the the neutral tions for </>_ and (J)n as one would expect. At finite tem- 
mode we derive the following contour ordered propagator peratures, Eq. (|B5() becomes 
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^l).!)' \'^') ^j 



In 






(B8) 



However, since the neutral mode is a non-propagating 
mode its propagator is independent of the temperature. 
Finally, in the limit of a; ^ 0^ and r ^ 0, Eq. HB8() 
and Eq. I|B7|1 becomes, 



G,,,v(0-,t) 



In 



sinh^ 






9 ^ 



G™^,(0-,i)= lim 



'Ujv^O+ 



.Itt 



^*-9X^,^'(0sgn(wArt) 

9 ^ 



(B9) 



-^-9X'),')'Wsgn(i) 



where we used 



In(if) = ln(i sgn(<)) + In \t\ = In |i| + i- sgn(i) (BIO) 

and lim,„j^^o+ sgn(wArt) = sgn(i). Also it is quite clear 
from the derivation that Gri^n'{x,t) = G,,',,,(— x, — t). 



APPENDIX C: MULTIPLE VERTEX FUNCTION 
CORRELATOR 

Consider the vertex operators V{q,r) = e'^i'l'i'^) , Here 
we calculate the correlation function of N vertex oper- 
ators (F(gi,ri) • • • F(gAr,rJv)) with respect to the free 
boson Lagrangian Eq. (|Aip 



(F(gi,fi)...F(gA,,nv)) 



P^e 



^So[4>] g»Ef=i9j0(f'j) 



Vcj) e 



-So[' 



(CI) 



Clearly this correlator will have the symmetry Eq. (|A2I 
only when 



N 



(C2) 



which is a charge neutrality condition whose important 
consequence will be shown below. Introducing a source 
term 






(C3) 



one can calculate the correlator using path integral as 
follows: 



fv(j)e 



iSo[<p] gi J d^ra(r)if>{r) 



V4) e 



iSo[0] 



(C4) 



p~^ f d^r f d^r' a{f)G{r—r^)a{f') 

5S,"fc = l9j9fcG(r,-rfc) 



= e 2 Z^3 

However the exponent can be separated into two parts 
since 

JV 

E IjIkGirj - fk) 

j,k=l 

N 

= E lMG{r, - r,) ~ G(0) + G(0)] 

3,k=l 

N N 

= iY, qjfG{0) + ^ qMG{rj - r,) - G(0)] 

N N 

= {Y,qjfG{Q) + 2 ^ q,qk[G{r, - fk) - GiO)]. 

j=l ]>k=l 

(C5) 

Notice though G(0) has IR divergence, as it was pointed 
out in Appendix combined with the charge neutrality 
condition Eq. IIC2p the correlation function becomes 

(y(gi,fi)...y(gw,rjv)) =e-^">.--i«^«''^(^-^-^"'') (C6) 

where the finite subtracted boson propagator G(r) = 
G{f) — G(0) naturally shows up. 



APPENDIX D: UNITARY KLEIN FACTORS 

Here we derive the unitary Klein factor algebraS used 
in Section IHII and their contour ordering. Anyonic ex- 
change statistics between qp's of different edges I and m 
defined in Eq. (|3.4|) requires 



t,/,t 



V'/V- 



I Vm 



= ±*8,/,t ,/,t 



^Ini^l- 



(Dl) 



However since boson fields commute between different 
edges, the statistical phase should come solely from the 
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Klein factor algebra. Hence if wc define real parameters 
aiyn by 



F,F — p-*"'"" F F: 



(D2) 



the magnitude of aim follows immediately from compar- 
ing Eq. IjDip with the following equation 



V^/^1„ = e— ^1„V/ 



(D3) 



to be 

ai„i = ±9 = -ami (D4) 

leaving the sign of aim still arbitrary. To consistently 
—--< < . , ^ <^ 




(a) 



(b) 



FIG. 10: Two configurations with the same topology: 
(a)Three edges are disconnected and each are described by 
If I with / = 0, 1, 2. The tunneling to edges I = 1,2 occurs at 
{xo = —a/2 = Xi , a;i =0) and {xq — a/2 = X2 , 2:2 = 0) respec- 
tively, (b) All three edges are connected with an open contour 
starting at 2; = which is described by a single chiral boson 
ifi. For the tunneling processes to edges I = 1,2, qp's hop out 
from X = —a/2 = Xi and x = a/2 = X2 respectively. We note 
in passing that since the choice of the origin is arbitrary and 
do not affect the end result, we made a convenient choice. 



J 



fix the correct sign, one must compare the system which 
is equivalent to Fig. ^| (a) with the same topology but 
where edges are connected with an open contour that 
brings all the edges into one as in Fig. llOf b) following 
Guyon et al.— . Klein factors are not necessary in the 
latter system because the commutation relations of tun- 
neling operators are enforced by the single chiral bosonic 
field 0. 

One way of determining Klein factors for case (a) is to 
compare the relative phases e"*"'" of Fig.^](a) with the 
phases e'^'^^'sniXi-x^) of Fig.[Tni(b) where Xi and Xm are 
now abscissas of the same edge measured from the point 
a; = (see the caption of Fig. ^| ). Another way is to 
look at the commutation relations of different tunneling 
operators Vi and to require that the commutator of the 
tunneling operators of the disconnected system give the 
same results as the connected one. 

These two approaches give the same results for the 
Klein factors and we will choose to detail the second ap- 
proach here. For the connected system of Fig. ^| (b), 
there is no need for the Klein factors since the tunnel- 
ing operators at different positions will commute due to 
chirality (as long as the tunneling paths do not cross). 
Hence we will set the constraint for a's by requiring Vi 
and Vz in the case of interest Fig. ^| (a) to commute. 
First, by comparing the exchange of ^(x)^{x') with that 
oiilj'^{x)'4){x') using the CBH formula, one finds that the 
definition of Eq. (|D2() implies 



FiF-^ = e^"'" 



F^^Fi. 



(D5) 



With Eqs. IJD2|I and IjDSp we can exchange Vi and V2 at 
equal time as the following 



Y^Yz = FQF~^FoF~^e^'^''^'^''^^^'^e^'^°^^°^^'^^e~'''''^^'^^^°'>e~^'^^^^^^°'> 

= FnF~^ FoF~^ 6^''°'°^'^°''^"'^°''^^ e^^ ^^^^^'^~^'^^ e^'^"^^^^ e^'^"^^'^^ e~'^''''^'^^ e~^'^'^'''^^ 



(D6) 



r 



where we used the fact that [ipi, tpm] = for I ^ m and straint: 



i<Po{xo=Xi} iipo{xo=X2) __ ^i9 sgn(Xi-X2) iipo{X2) i<poiXi) 

(D7) 
The second equality followed from qp exchange within 
edge given by Eq. H2.6|l and series of Klein factor ex- 
change as the following: 

FnF^^FoF^^ = F^F^^ F^^ Foe'""^ 

= FoF^'F^'Foe" 






(D8) 



From Eq. ljD6p . requiring Vi and V2 to commute as they 
would in a single connected edge, results in the con- 



ao2 + "21 + "10 + & sgn(Xi - Ar2) = 0, 



(D9) 



which allows us to define a^z — "21 = ^"10 — ^ since 
Ar2 — ATi =a> 0(See the caption of Fig. I1U|) . Further we 
can conclude 



FqFi F0F2 — F0F2 FqF^ e' . 



(DIO) 



Now we determine the contour ordered Klein factor 
"propagator" {TKiFoF^^Y{t^)iFoF^^Y' it'^'))o based 
on Eq. I|D10() . We first note contour ordering means that 
(FoFf i)'(t'') and {FoF-'^y' {t"^' ) need to be exchanged 
whenever t"' >k t''', i.e. 
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{TK{FoF-'r{f^){FoFi'r'{t'^'))o 



{{FoF,-'Y{t^){FoF^'Y'{t"^')) for t^ >k i'"' 
{iFoF^'y'it'^'')iFoF^'yin) for t'^'>K i" 
I 



(Dll) 



From Eq. IJD10() , it is clear that a consistent way of defin- 
ing the contour ordered propgator is to assign a relative 
phase factor of e-'^^^snlXi-Xs) between t"''>K f^ case (in 
which the Klein factors need to be exchanged for contour 
ordering) and f >x i"* case(in which no exchange is 



neccesary). Further in order for {TKViit'^)V2it''^' )) = 
{TKV2{t"'')Vi{t'i)) to hold, which would be the case for 
the connected edge and which connects us back to the 
basis for the derivation of Eq. IJD10|) , we require 



{TKiFoF,-'nf')iFoF,-'y'it'^^')) = {TK{FoF,-'r'it'^^'){FoF,-'Y{t^)). 

I 



(D12) 



This last requirement can be satisfied by assigning half of Eq. IJD11(I and we arrive at the following form: 
of the relative phase factor e~'*'^sn(-^i--f2) ^q each case 



{TK{FoF^y{t'^){FoF^y{t"^'))o = e'"'f-s"(-^i-^2)x..v (*-*') ^ e-«"'|x„,v (*-*') (D13) 

I 



where we used the definition of Xv.n'i^ ~ ^') introduced Eq. IBGI 
for the contour ordered propagator of chiral bosons in 
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